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山东 省 植被 NPP 时 空 分 布 特征 及 驱动 因素 分 析 
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摘要 : 为 了 解 区 域 生 态 系统 植被 生产 力 的 分 布 规律 及 其 变化 对 


的 年 NPP 数据 ， 采 用 一 元 线性 


区 域 生 态 调控 的 作用 ， 本 文 基于 MOD17A3 
省 柑 


回归 分 析 法 和 Pearson 相关 系数 法 等 方法 ， 分 析 了 2000 一 2014 年 山东 


被 VPP 的 时 空格 


局 变化 特征 及 其 气候 、 人 为 等 


植被 NPP 总 体 呈 上 升 趋势 ， 


间 分 布 具有 明显 
布 于 山东 省 西部 
占 山东 省 植被 总 


的 分 异 特征 ， 其 
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积 的 92.63%。 


A 响 。 


年 均值 为 442gC.m*.a， 
中 年 均 NPP 300 ~ 400gC.m 
的 绝 大 部 分 区 域 。(2) 山东 省 植被 NPP 年 增加 速率 平均 值 为 3.29gC.m? 
其 中 ， 极 显著 增加 区 域 占 山东 省 总 


es 
年 总 量 


结果 表明 : 
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“a 的 


平均 值 为 63.16Tg C。 山 东 


再 积 最 大 ， 


积 的 16.56%， 主 要 


青岛 市 、 


日 照 市 、 


询 沂 市 


\ 潍 坊 市 


\ 东 营 市 和 滨州 市 ; 显著 减少 和 极 显 著 


(1) 2000 一 2014 年 ， 山 东 省 
省 植被 NPP 空 
占 总 面积 的 52.28%， 分 
. a.， 增 加 区 域 
分 布 在 


因 台 市 、 


减少 的 区 域 占 山 


东 省 总 面积 的 1.33%， 


js 


符 与 


南部 和 北部 黄河 三 角 洲 地 区 。(3) 植 被 NPP 受气 候 变化 和 人 类 活动 等 因素 的 共同 影响 ， 
气候 因素 中 降雨 的 影响 最 大 ， 城 市 化 随 着 人 类 活动 的 频繁 逐渐 成 为 植被 VPP 变化 的 重要 驱动 力 。 

关键 词 : NPP， 年 际 变化 ， 驱 动因 子 ， 山 东 省 ， 土 地 利用 类 型 
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Temporal-spatial variation characteristics and driving factors of 


分 布 在 山东 省 东 


[mi 


vegetation NPP in Shandong Province 


LUO Yan, ZHANG Songlin 
(Northwest Normal University, Lanzhou 730070, China) 
Abstract: In order to understand the distribution law of regional ecosystem vegetation productivity and its effect 
on regional ecological regulation, based on the annual net primary productivity (NPP) data of MOD17A3, the 
method of univariate linear regression analysis and Pearson correlation coefficient, etc., were used to analyze 
spatial and temporal patterns of vegetation NPP ecosystem in Shandong Province during the past 15 years from 
2000 to 2014 and climate and human factors. The results showed that: (1) the average NPP of vegetation in 
Shandong Province from 2000 to 2014 was 442g C ， m”. al and the total average annual NPP was 63.16TgC. 
The vegetation NPP of Shandong Province showed obvious spatial differentiation characteristics, among which, 
the vegetation area of the average NPP 300 ~ 400 gC.: m’”. al! was the largest. (2) the average annual growth 


rate of NPP in Shandong province was 3.29 g C ma accounting for 92.63% of the total vegetation area in 
Shandong province. Among them, the area of extremely significant increase accounted for 16.56% of the total 


area of Shandong province, mainly distributed in Yantai, Qingdao, Rizhao, Linyi, Weifang, Dongying and 
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Binzhou cities; The regions with significant decrease and extremely significant decrease accounted for 1.33% of 


the total area of Shandong province, and were scattered in the southeast of Shandong province and the northern 


Yellow River delta region. (3) the vegetation NPP by the combined effect of factors such as climate change and 


human activities. The influence of rainfall on climate factors is biggest. Urbanization gradually becomes an 


important driving force for vegetation NPP changes with the frequent human activities. 


Key words: net primary productivity, inter-annual variability, driving factors, Shandong Province, land use type 


植被 净 初 级 生产 力 (NPP，Net Primary Productivity) 指 绿色 植物 在 单位 时 间 
机 干 物质 量 (Field et al.，1998), 体现 陆 
的 能 力 ( 张 杰 和 潘 晓 玲 ，2010)。 植 被 NPP 在 全 球 的 碳 平衡 中 具有 重要 的 作 月 
科学 研究 的 焦点 之 一 《Cramer et al.，1999)。 近 和 


态 系 统 当时 的 状况 ， 


面积 上 所 积累 的 有 


飞 表 绿 色 植 物 通 


究 。 李 登科 等 (2018) 基于 MOD17A3 数据 对 


变化 的 关系 进行 研究 , 发 现 海拔 高 度 不 是 主要 的 影响 


着 快速 
致 了 植被 及 水 土 资源 的 流失 , 研究 NPP 的 时 3 
省 生态 环境 起 着 重要 作 月 
定 和 碳 循环 的 调节 有 具 有 重要 芯 
目前 针对 山东 省 植被 NPP pp 的 相关 和 
空 分 布 特征 及 其 气候 、 人 为 等 
源 的 有 效 合 理 利 用 提供 科学 依据 。 


山东 


NPP 的 影响 。 
植被 NPP 的 时 


的 区 域 VPP 值 越 高 。 
山东 省 植被 种 类 繁多 ， 并 且 作 为 中 国 的 农业 大 省 ， 
束 的 经 济 发 展 ， 人 类 活动 的 干扰 逐渐 加 强 ， 使 得 土地 收益 面积 
变化 可 以 有 效 地 反映 其 生态 
日 。 SO 与 气候 变化 以 及 相关 影响 因 


区 变化 趋势 不 明显 ， 呈 现 显 
型 反 演 了 全 国 植被 WPP， 认 为 近 十 各 
岂 域 性 明显 ， 且 降雨 量 和 温度 为 影响 植被 NPP 的 主要 
类 型 和 气候 变化 的 影响 表现 出 较 大 的 时 空 和 异 质 性 ， 
进一步 研究 。 陈 瑛 瑛 等 (2015) 分 析 了 我 国 


著 增 加 趋势 和 显著 减 
FE 来 中 国 植被 NPP 的 空 


过 光合 作用 固定 大 气 CO， 
日 ， 对 气候 变化 的 响应 研究 是 


FE 来 ， 多 位 学 者 对 我 国 植被 NPP 与 气候 因子 的 关系 进行 了 
国 2000 一 2015 年 陆地 植被 NPP 进行 研究 发 现 ， 大 部 分 
少 趋势 的 面积 相对 较 少 。 仲 晓 春 等 〈2016) 利用 CASA 


间 分 布 总 体 北 低 南 总 
因素 。 然 而 ， 在 全 国 尺度 上 ， 植 被 NPP 受 不 同 植被 


因此 植被 NPP 在 区 域 尺度 上 的 


西北 向 东南 增加 ， 


时 空 变化 的 驱动 机 制 需要 
净 初 级 生产 力 ， 认 为 大 部 分 地 区 植被 NPP 与 
温度 呈正 相关 、 与 降雨 呈 负 相关 ， 受 气候 因素 影响 大 。 王 志 红 等 (2018) 对 贵州 植被 NPP 分 布 特点 与 地 形 
因子 ， 地形 起 伏 度 与 NPP 显著 相关 ， 且 地 形 起 伏 越 大 


、 玉 米 和 花生 等 作物 的 主要 生长 区 。 然 而 随 


发生 改 变 ， 土 地 利 ) 


] 类 型 发 生 转 移 ， 导 


究 较 少 ,， 本 文 利 用 MODIS 数 和 
因素 对 其 的 影响 ， 为 山东 省 的 植被 保护 、 碳 循环 评估 和 自然 资 


1 研究 区 概况 


地 


山东 省 地 理 坐 标 114*47' EE 一 122°43' EE， 34°22' N 一 38°%27'N, 位 于 我 
黄海 和 渤海 ， 陆 地 总 面积 15.8 2 km’, 


10 006 万 人 【2017 年 )， 


也 ， 北 部 和 西部 为 平坦 的 黄 ; 


冲积 平原 ， 东 部 与 南部 为 起 伏 的 丘陵 区 。 


0 
1 


系统 的 动态 和 


[健康 状况 ,对 评价 
素 的 研究 对 生态 系统 的 稳 
究 了 中 国 北方 因 土 地 利用 和 履 盖 方式 的 改变 对 于 
中 和 GIS 相关 技术 分 析 山 东 


国 东 部 沿海 地 区 和 黄河 下 游 ， 东 临 
{£17 个 地 级 市 。 中 部 为 隆起 的 山 


自 较 大 ， 由 季 分 明 ， 雨 量 集 


2 数据 来 源 与 处 理 


本 文 使 用 的 2000 一 2014 年 的 MOD17A3 NPP 数 
(http://Ipdaac.usgs.gov)， 空 间 分 


立 起 全 球 NPP 估算 模型 ， 


为 期 15 年 的 MODIS 数据 在 MRT 工具 


中 ， 夏 季 雨 热 同 期 ， 冬 季 了 晴朗 干燥 ， 平 均 气 温 


00 mm， 从 东南 沿海 向 西北 内 陆 递 减 。 


该 数据 是 也 
人 NPP。 目 前 ， 该 数据 在 全 


i 


的 时 空 特征 和 碳 循环 研究 


I 


oi 


本 


合 1 000 m 的 DEM 数 


到 不 同 地 


暧 温带 季风 气候 ， 受 海洋 季风 的 影 
11 ~ 14 '‘C; 年 均 降 水 量 500 ~ 1 


居 来 自 美国 NASA EOS/MODIS 网 站 
上 用 BIOME 一 BGC 模型 与 光 能 利用 模型 建 
x 植被 净 初 级 生产 力 


2009; Rasmus et al.，2006; 位 贺 杰 等 ，2014)。 本 文 将 
\ 行 处 理 ， 投影 转换 为 Albers Equal Area。 

究 使 用 的 气象 数据 来 自 中 A (http://cdc.cma.gov.cn) 的 山东 省 境内 及 周 
边 的 33 个 气象 站 点 提供 的 年 均 气 温和 年 降水 量 数据 。 


中 和 33 个 气象 站 点 的 经 续 


度 信息 ， 利 用 Aunsplin 气象 数据 插值 软件 对 气象 数据 进行 插值 处 理 ， 获 取 的 栅 格 分 辩 率 为 1 kmx1l km。 
Aunsplin 软件 对 气象 数据 插值 去 除了 高 程 对 温度 和 降水 的 影响 ,得 到 的 气象 栅 格 数据 更 加 可 靠 ( 程 春晓 等 ， 
2014)。 

土地 利用 类 型 数据 为 2000 年 、2005 年 和 2010 年 山东 省 遥感 监测 数据 , 来 源 于 中 国 科 学 院 资源 环境 科 


学 数据 中 心 , 分 辩 率 为 1 km。 将 土地 利用 类 型 数据 归并 分 类 为 林地 、 耕 地 、 水 域 、 草 地 、 城 镇 与 未 利用 地 。 
2.1 NPP 变化 趋势 分 析 
山东 省 2000 一 2014 年 的 平均 NPP 空间 分 布 图 ， 计 算式 为 : 
NPP = Ler (1) 
式 中 ，NPP 为 山东 省 2000 一 2014 年 累计 15 年 的 NPP 平均 值 ; n 为 研究 时 间 段 的 年 数 ， 取 值 15; 为 
年 序号 ， 即 为 2000，2001，...，2014 年 ;NPP 为 第 年 的 山东 省 NPP 值 。 


基于 像 元 的 NPP 线性 倾向 采用 一 元 线性 回归 分 析 法 ， 消 除了 异常 因素 对 NPP 变化 趋势 的 影响 〈 潘 竟 
虎 和 李 真 ，2015)。 本 文 对 NPP 年 际 变 化 趋势 的 计算 公式 为 : 


NXP Pid i 

式 中 ，Slope 为 线性 倾向 值 ， 单 位 为 gC . m?. a!， 代 表 了 NPP 的 变化 倾向 程度 。 当 Slope > 0，NPP 
随 着 时 间 增 加 呈 上 升 趋势 ， 当 Slope = 0，NPP 随 着 时 间 增 加 不 变 ; Slope < 0，NPP 随 着 时 间 增 加 呈 下 降 
趋势 。n 为 研究 时 间 段 的 年 数 ， 取 值 15; i 为 年 序号 ， 即 为 2000，2001，...，2014 年 ， NPP 为 第 年 的 NPP 
值 。 


Slope = (2) 


可 归 模 型 的 显著 性 大 小 采用 于 检验 ， 计 算式 为 : 
_ TEPPi NPP) 
Fie1(NPPi-NPP)? 
式 中 ，NPP; 为 第 年 的 NPP 值 ，NPP 为 累计 15 年 NPP 的 平均 值 ，NPP 为 第 :年 的 NPP 回归 值 。 
根据 F 检验 查询 结果 ， 将 山东 省 NPP 的 变化 趋势 划分 为 极 显著 减少 、 显 著 减 少 、 变 化 不 明显 (不 显 

著 增 加 或 者 减少 )、 显 著 增 加 和 极 显 著 增 加 五 个 类 型 。 

2.2 相关 性 分 析 
基于 像 元 的 植被 NPP 与 气候 因子 年 均 气温 和 年 降水 量 的 相关 性 分 析 采 用 Pearson 相关 系数 法 。 之 后 对 
NPP 与 气温 或 降水 的 相关 系数 进一步 计算 ， 得 到 NPP 与 气温 和 降水 的 偏 相关 系数 。 相 关系 数 以 及 偏 相 关 
系数 的 计算 方法 如 式 (4) 和 式 (5) 所 示 : 


(n—2) (3) 


Bia[CXi—X) (yi—y)] 
EDxT 0 D 
式 中 ，Ry 为 两 变量 NPP 与 气温 或 降水 的 相关 系数 。 
Ryz — Rxy-RxzRyz __ (5) 
|a- RDxG-R3a) 
式 中 ，Rxyz 是 气温 或 降水 〈z) 固定 后 另 一 气候 因子 与 NPP 的 偏 相关 系数 。 
偏 相 关 性 的 显著 性 检验 采用 t 检 验 ， 计 算式 为 : 
t=Vn -m1 (6) 


JR2 yz 


(4) 


式 中 ，n 为 年 份 ，m 是 自 变 量 个 数 。 

各 个 要 素 之 间 的 相互 影响 较为 复杂 ， 仪 研究 各 要 素 间 的 相关 性 和 偏 相 关 性 不 能 反映 出 他 们 之 间 的 关系 ， 
复 相 关 性 可 以 进一步 反映 各 要 素 间 的 综合 影响 ， 因 此 基于 相关 系数 和 偏 相 关系 数 计 算 复 相关 系数 ， 公 式 如 
下 : 


Rx,yz = /1 一 (1— R#y)(1 一 R22y) 


式 中 ，R%yz 为 因 变 量 x 与 自 变 量 人 z 的 复 相关 系数 ，Re, 为 俩 相关 系数 ，Rxy 为 x 与 y 的 相关 系数 。 
复 相关 系数 的 显著 性 检验 (下 检验 ) 公式 为 : 


Rxyz X Rxyz 有 n=k=1 
1 — Ryz X Rx,yz k 


F= 


式 中 ，n 为 年 份 ,kk 是 自 变量 个 数 。 
2.3 气候 因子 与 非 气候 因子 的 驱动 力 分 类 
参考 国内 学 者 的 研究 〈 高 志 强 等 ，2004; 何 奕 营 等 , 2018) 和 植被 覆盖 变化 驱动 分 区 原则 (陈云 浩 等 ， 
2001)， 并 结合 研究 区 的 情况 对 其 进行 适当 调整 ， 把 影响 NPP 变化 的 驱动 力 分 为 五 类 〈 表 1)。 
表 1 NPP 变化 驱动 力 分 区 


Table 1 Classification criteria for NPP change driving force 


变化 驱动 因子 分 区 准则 Rules 
Changes driving factors RR RB 

气温 降水 强 驱 动 [T+P]* >t00 (>t00 F>Foos 

气候 因子 气温 降水 弱 驱 动 [T+P] lwo lstoo F> Foos 
Climate factors 气温 为 主要 驱动 T ld>toor F> Foos 
降水 为 主要 驱动 P lf>t00 F> Foos 

非 气候 因子 E 加 

非 气候 因子 驱动 NC F<Foos 


Non-climate factors 
注 : R1 为 忽略 温度 影响 的 偏 相关 性 显著 性 检验 ，R2 为 忽略 降水 影响 的 偏 相 关 性 显著 性 检验 ，R3 为 复 相 关 的 显著 性 检验 


Note: RI1 is a partial correlation significance test that ignores temperature effects; R2 is a partial correlation significance test that 


ignores precipitation effects; R3 is a significant test of complex correlation 


驱动 力 分 析 的 统计 检验 方法 如 下 : 

对 于 R1 和 R2 偏 相 关 性 的 显著 性 检验 采用 t 检 验 , 当 显 著 性 水 平 设 定 为 0.01 时 , t 取 值 为 2.681, 当 上 > 
2.681 时 ， 显 示 显 著 。 对 于 R3 复 相 关 性 的 显著 性 检验 采用 下 检验 ， 当 显著 性 水 平 设 定 为 0.05 时 , FF 取 值 为 
3.885， 当 五 > 3.885 时 ， 显 示 显 著 。 

非 气候 因子 中 十 地 及 NPP 转移 矩阵 的 获得 方法 为 : 

将 2000、2005 和 2010 年 三 期 土地 利用 数据 在 GIS 软件 中 分 类 得 到 林地 、 耕 地 、 水 域 、 草 地 、 城 镇 与 
未 利用 地 六 种 土地 利用 类 型 ， 并 赋予 不 同 的 值 。 把 2000 年 和 2005 年 两 期 土地 利用 数据 相 减 ， 得 到 2000 
一 2005 年 不 同类 型 土地 向 男 外 五 种 土地 类 型 转移 的 栅 格 图 ， 其 中 每 一 种 转移 方法 都 有 特定 的 值 ，2005 一 
2010 年 的 方法 同上 ， 最 后 进行 统计 ， 得 到 土地 转移 和 矩阵。 利用 得 到 的 土地 利用 类 型 转移 数据 分 别提 取出 
2000、2005 和 2010 年 的 NPP 数据 ， 进 行 计 算 ， 得 到 2000 一 2005 年 和 2005 一 2010 年 各 土地 类 型 转移 的 
NPP 量 ， 统 计 得 到 NPP 转移 矩阵 。 


3 结果 与 分 析 


3.1 植被 VPP 的 空间 分 布 特征 

山东 省 2000 一 2014 年 植被 的 平均 WPP 为 442 g C . m2 .a!。 由 图 1 可 以 看 出 ， 山 东 省 植被 NPP 平均 
值 为 0~200gC.m?.a! 的 区 域 很 少 ， 故 合并 为 一 组 ，200 ~ 400 gC.m?*. a 的 区 域 分 布 较 集中 ， 以 50 g 
C. m”. a 为 级 差 可 有 效 区 分 其 空间 分 布 ，400 gC. m? . a 以 上 的 分 布 较 离散 ， 以 100 gC.m?. a- 为 级 
差 即 可 有 效 区 分 其 空间 分 布 。 山 东 省 植被 NPP 年 均值 空间 分 布 呈 现 明显 的 分 异 特征 ， 其 中 ， 平 均 NPP 低 
于 300 gC. m* .al 的 区 域 集中 分 布 于 济南 市 、 淄 博 市 、 莱 芜 市 、 泰 安 市 和 济宁 市 ， 在 德州 市 西南 部 、 滨 
州 市 和 东营 市 北部 地 区 也 有 分 布 ， 占 总 面积 的 11.72%; NPP 平均 值 300 ~ 400 g C . m?. a! 的 分 布 于 山东 
省 西部 的 绝 大 部 份 区 域 ， 占 总 面积 的 一 半 以 上 (52.28%); NPP 年 均值 400 ~ 500 gC . m*. a 的 多 数 分 布 


并 


| 


201903.00206v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


于 潍坊 市 、 日 照 市 和 临沂 市 ， 少 数 分 布 于 菏泽 市 、 束 庄市 南部 、 以 及 烟台 市 和 青岛 市 的 西部 ， 占 总 面积 的 
13.89%; NPP 平均 值 大 于 500 gC.:m*…al 的 区 域 位 于 山东 省 东北 部 沿海 区 域 ， 占 比 为 22.11%， 最 大 值 出 现 
在 威海 市 。 鲁 中 南山 地 丘陵 区 NPP 年 均值 为 250 ~ 300 gC . m”. al1， 主 要 是 由 于 地 形 抬升 、 降 水 减少 ( 张 
可 欣 和 孙 成 武 ，2009) 和 和 气温 降低 所 致 ， 东 部 沿海 大 多 区 域 的 植被 NPP 年 均值 为 500 ~ 700 gC. m2? . an 
与 其 海拔 较 低 、 气 温 较 高 和 降水 相对 丰沛 有 关 。 


_NPP 平均 值 
NPP average 
[gC'm™a] 

国 400-500 
200 ” 国 500~600 
00~250 国 600-700 
50-300 国 :700 


1 2000 一 2014 年 山东 省 NPP 分 布 
Fig.1 Distribution of NPP in Shandong Province from 2000 to 2014 

3.2 植被 VPP 的 时 间 分 布 特征 

2000 一 2014 年 ， 山 东 省 植被 NPP 总 量 年 平均 值 为 63.16 Tg C， 随 着 时 间 延 续 有 升 有 降 ，2000 一 2004 
年 ， NPP 持续 增加 (图 2)， 之 后 ， 相 对 稳定 略 有 升 高 。2004 年 NPP 总 量 最 高 ， 为 70.59 T g C， 其 次 是 
2008 年 ， 为 69.74 Tg C， 两 者 的 NPP 值 分 别 高 于 平均 值 的 11.76% 和 10.42%; 2000 年 总 量 最 低 ， 为 53.34 
Tg C， 低 于 平均 值 15.55%; 最 低 值 仅 是 最 高 值 的 0.76 倍 ， 表 明山 东 省 植被 的 NPP 值 的 年 际 变 化 较 大 、 不 
太 稳 定 。2000 一 2002 年 , 植被 NPP 值 低 于 15 年 的 平均 值 , 主要 是 由 于 这 些 年 山东 省 连续 干旱 ( 张 胜 平等 ， 
2004), 使 部 分 区 域 的 植被 生产 力 降低 〈 杨 晓 霞 等 ,2003)。2005 年 山东 省 植被 NPP 又 下 降 到 61.79TgC， 
可 能 是 由 于 年 初 暴风 雪 天 气 影响 了 植被 的 生长 〈 崔 林 丽 等 ，2016)。 随 后 几 年 ， 植 被 生产 水 平 逐渐 恢复 ， 
NPP 总 量 总 体 略 微 升 高 。2009 一 2011 年 植被 NPP 总 量 较 低 ， 可 能 与 当时 气温 与 降水 都 减少 有 关 。 如 图 3 
所 示 ，2009 一 2011 年 研究 区 平均 气温 和 降水 量 都 呈现 减少 趋势 。 
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图 2 山东 省 2000 一 2014 年 VPP 总 量 
Fig.2 Total NPP in Shandong Province from 2000 to 2014 
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注 : A. 平均 气温 ，B. 累积 降水 量 。 
Note: A. Average temperature; B. Cumulative precipitation. 
图 3 2000 一 2014 年 山东 省 年 均 气 温和 年 降水 量 的 年 际 变化 


Fig.3 Interannual variation of annual average temperature and annual precipitation in Shandong Province from 


2000 to 2014 


2000 一 2014 年 ， 山 东 省 植被 NPP 的 Slope 平均 值 为 3.29 gC.m?*.a!。 图 4 表明 ， 山东 省 92.36% 的 
植被 NPP 呈 增 加 趋势 ，7.64% 的 植被 NPP 为 减 小 趋势 。 尺 检验 表明 ， 山 东 省 62.95% 的 植被 NPP 年 际 变化 
趋势 不 明显 ; NPP 年 际 变化 为 正 值 区 域 〈 增 加 趋势 ) 的 平均 为 3.93 gC. m*.a!; 极 显著 增加 区 域 占 山东 


省 植被 面积 的 的 16.56%， 主 要 分 布 在 烟 


台 市 、 青 岛 市 、 潍 坊 市 、 日 照 市 、 临 沂 市 、 东 营 市 和 滨州 市 ， 显著 


增加 区 域 占 山东 省 植被 总 面积 的 19.16%， 大 部 分 位 于 极 显著 增加 区 域 的 周边 ， 少 部 分 在 泰安 市 和 德州 市 ; 


NPP 呈现 减少 趋势 的 Slop 


零星 分 布 在 山东 省 东南 部 条 


e 平均 值 为 -2.80 gC . m”. a ， 其 中 显著 减少 和 极 显著 减少 的 部 分 仅 有 1.33%， 
[北部 黄河 三 角 洲 地 区 《图 4)。 


山东 省 不 同 土地 覆盖 类 型 的 年 际 增长 速率 各 不 相同 , 植被 NPP 在 草地 中 以 4.69 gC. m”. a 的 平均 速 
率 增长 , 且 增 长 趋势 最 明显 , 其 中 呈现 极 显著 增加 趋势 和 极 显 著 减 少 趋势 的 面积 分 别 占 草地 总 面积 的 26.43% 
和 0.37%; 林地 的 增长 速率 为 4.54 gC . m”. a', 仅 在 草地 之 下 , 呈 极 显著 增 和 极 显 著 减 的 面积 分 别 为 20.38% 


和 0.25%。 耕 地 、 城 镇 和 水 域 的 平均 增长 速率 分 别 为 3.39、2.37 和 2.60g C . m” . al1， 相 应 的 极 显著 增加 


面积 15.61%、12.17% 和 14.77%, 极 显著 减少 面积 0.53%、2.54% 和 1.11% 。 而 耕地 、 城 镇 和 水 域 变 化 不 明显 
的 区 域 与 草地 与 林地 相 比 相对 较 多 ， 分 别 为 64.59%、68.91% 和 65.45%， 可 见 ， 山 东 省 增长 速率 较 明 显 的 
草地 和 林地 ， 其 极 显著 增加 区 域 的 面积 占 比 较 大 ， 而 显示 的 极 显著 减少 的 面积 相对 较 少 ， 无 明显 变化 的 区 


域 面积 占 比 小 。 
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图 4 2000 一 2014 年 山东 省 NPP 年 际 变化 趋势 和 显著 性 分 布 
Fig.4 Interannual variation trend and significance distribution of NPP in Shandong Province during 2000 to 
2014 


3.3 植被 NPP 变化 的 驱动 因素 
3.3.1 气象 因子 对 植被 WPP 的 影响 

植被 NPP 受气 候 变 化 的 影响 很 大 , 分 析 气 候 变化 特征 有 助 于 进一步 探讨 植被 NPP 的 影响 因素 (Gao et 
al.，2013)。 如 图 5， 山 东 省 植被 NPP 值 与 年 均 气 温 的 平均 相关 系数 为 0.09， 表 现 为 很 弱 的 正 相 关 关 系 。 
正 相关 区 域 面积 占 植被 总 面积 的 70.58%。 负 相关 区 域 面积 占 总 面积 的 29.40%， 了 检验 发 现 ， 显 著 相 关 
P<0.01) 的 区 域 仅 有 0.16%, 位 于 日 照 市 的 沿海 地 区 。 植被 NPP 年 均值 与 年 降水 量 平均 相关 系数 为 0.28， 
较 气 温 与 NPP 的 相关 性 略 强 。 正 、 负 相关 区 面积 分 别 为 91.49% 和 8.50% 。 了 检验 可 知 ， 有 1.49% 的 区 域 通 
过 了 P< 0.01 的 显著 性 检验 ， 主 要 分 布 在 鲁 中 山地 ， 零 散 分 布 在 山东 省 黄河 三 角 洲 等 地 。 相 比 之 下 ， 山 东 
省 植被 NPP 与 年 均 降 水 量 的 相关 性 更 强 且 更 显著 。 


后 浊 相 关系 区 ? 所 降水 相 立 

靳 。“ 温 相关 系数 WY 降水 相关 系数 

3 Temperalure ) 2 Precipitation 
correlation coefficient 人 correlation coellicient 


高 High:0.78 


4 
| 让 l £ 作 ) -二 
低 Low:-0.58 a y 


低 Low:-0.63 
口 省 界 


Provincial boundary 


图 5 植被 NPP 与 气候 因素 相关 系数 分 布 
Fig.3 Distribution of correlation coefficient between vegetation NPP and climatic factors 
不 考虑 降水 的 影响 ,植被 NPP 与 温度 年 均值 的 平均 偏 相关 系数 为 0.24， 显著 相关 区 域 占 1.00%， 比 考 
虑 降水 影响 时 的 面积 占 比 增 大 了 0.83%; 忽略 温度 的 影响 ,植被 NPP 与 年 降水 量 平均 值 的 平均 偏 相关 系数 
为 0.35， 而 显著 正 相 关 区 域 占 6.54%， 比 考虑 年 均 气 温 影 响 时 的 面积 占 比 增加 了 5.05%， 不 显著 正 相 关 性 
区 域 减少 了 3.62%; 这 说 明山 东 省 植被 NPP 受 水 热 条 件 的 共同 作用 ,其 中 降水 较 气温 的 影响 更 重要 (图 6)。 
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图 6 植被 WPP 与 气候 因素 的 偏 相关 系数 分 布 
Fig.6 Distribution of partial correlation coefficient between vegetation NPP and climatic factors 
图 7 可 见 ， 复 相关 关系 在 0.003 ~ 0.881 之 间 。 植 被 NPP 与 气候 因子 复 相关 性 较 强 的 区 域 分 布 在 菏泽 
市 、 济 宁 市 和 现 庄 市 一 带 ， 复 相关 性 较 弱 的 区 域 位 于 烟台 市 、 青 岛 市 和 潍坊 市 。 植 被 NPP 与 气候 因子 的 地 
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区 差异 性 可 能 与 海拔 有 关 ， 青 岛 市 、 潍 坊 市 和 日 照 市 海拔 较 低 ， 而 周边 海拔 升 高 ， 低 海拔 地 区 气温 较 低 ， 
使 得 植被 NPP 与 气候 因子 的 复 相关 性 较 弱 。 
植被 NPP 受气 温 、 降 水 等 气候 因子 的 变化 与 人 类 活动 因素 的 影响 ， 降 水 对 其 的 影响 最 大 ， 而 人 类 活动 


中 ， 土 地 利用 类 型 的 转变 改变 了 生态 系统 的 类 型 和 结构 ，i 


而 改变 了 植被 生产 力 。 


由 图 7 驱动 类 型 图 可 见 , 山东 省 面积 的 0.56% 以 气温 降水 为 强 驱 动力 ,主要 位 于 在 菏泽 的 西部 和 枣庄 的 
北部 ， 以 气温 降水 弱 驱 动 的 区 域 大 多 分 布 在 山东 省 南部 几 市 ， 分 布 较为 分 散 ， 占 总 面积 的 1.76%; 气温 驱 
动 型 占 研究 区 面积 的 0.97%， 主 要 集中 在 菏泽 市 和 枣庄 市 ， 降 水 驱动 型 占 5.91%， 大 致 分 布 于 北部 的 滨州 


市 和 东营 市 ， 以 及 南部 的 泰安 市 、 菏 泽 市 、 


济宁 市 与 枣庄 市 的 交界 处 ; 非 气 候 因 子 驱 动 型 的 面积 达到 了 


90.79%， 在 山东 省 各 市 均 有 分 布 。 由 此 可 以 看 出 ， 非 气候 因子 对 山东 省 大 多 区 域内 的 植被 NPP 变化 影响 
较 大 ， 是 主要 的 驱动 因素 。 
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图 7 植被 NPP 与 气候 因素 的 复 相 关系 数 分 布 


Fig.7 Distribution of complex correlation coefficient between vegetation NPP and climatic factors 


3.3.2 土地 利用 变化 对 植被 NPP 的 影响 


如 表 2 所 示 : 2000 一 2005 年 有 3 106 km 的 土地 利用 类 型 发 生 了 转移 , 占 土地 利用 总 面积 的 2%， 其 中 


耕地 转移 为 最 主要 的 转移 类 型 , 转移 面积 
年 山东 省 共有 1 426 km 的 土地 向 其 他 土地 类 型 发 生 了 转移 ， 占 总 面积 的 0.92%, 耕地 转移 面积 占 总 转移 面 
积 的 70.48%， 主要 向 城镇 月 


总 转移 面积 的 51.09%, 主要 转移 为 城镇 和 水 域 用 地 。2005 一 2010 


地 转移 。 山 东 省 近 十 年 主要 在 耕地 与 城镇 、 未 利用 地 与 城镇 之 间 土 地 类 型 发 生 


转移 ， 主 要 是 以 耕地 、 林 地 、 草 地 和 未 利用 地 转 出 ， 以 及 城镇 和 水 域 用 地 转 入 为 主 的 土地 类 型 转移 。 
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土地 利用 类 型 
Land use type 
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Grassland 
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Water area 
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Urban area 


未 利用 土地 


表 2 土地 面积 转移 矩阵 


Table 2 Transfer matrix of land area (km’) 


耕地 林地 草地 水 域 城镇 未 利用 土地 
Farmland Forest Grassland Water area Urban area Unused land 
101 896 13 0 313 1 248 13 
4 9736 0 12 27 0 
481 17 13 163 109 81 26 
58 3 49 5641 146 14 
18 0 2 11 19 953 1 
185 1 3 44 227 1 834 
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向 较 高 的 的 转换 ， 使 得 植被 NPP 增加 (3 
东 省 植被 NPP 转移 的 情况 如 表 3 所 示 : 
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Table 3 NPP transfer matrix (x10*4 TgC) 


耕地 


Farmland 


24.07 


13.90 


114.42 


23.66 


表 3 NPP 转移 矩阵 

林地 草地 
Forest Grassland 
7.53 0 

0 
1.19 a 

0 0 

0 1.07 

0 0.16 
4.36 0 

和 0.43 

0 水 
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日 类 型 产生 的 NPP 各 有 不 同 。 不 同 土地 利用 类 型 间 植 被 NPP 的 转移 主要 表现 为 具有 
也 履 盖 类 型 向 较 低 NPP 的 的 土地 有 覆 六 类 型 的 转换 , 使 得 植被 NPP 损失 ; 或 是 植被 NPP 
E 芳 等 ，2018)。2000 一 2010 年 ， 不 同 土地 利 


未 利用 土地 
Unused land 


6.31 


0.16 


未 利用 土地 
Unused land 
注 :“*” 代 表 该 方式 下 不 转换 。 


Note: “*” means no conversion in this mode. 


2000 一 2005 年 ， 草 地 和 未 利用 地 转移 为 耕地 使 得 NPP 增加 量 最 大 ， 为 386.7x10“T g C; 耕地 转移 为 
城镇 用 地 使 植被 NPP 减少 156.69x10* Tg C; 各 土地 覆盖 类 型 植被 NPP 向 草地 转移 的 WPP 量 最 少 ， 仅 有 
1.23x10*TgC， 其 次 为 林地 ， 为 8.72x104TgC， 耕 地 最 多 ， 有 426.54x104Tg C。 

2005 一 2010 年 ， 城 镇 向 耕地 和 水 域 用 地 的 转换 ， 使 得 植被 NPP 增加 量 最 大 ， 为 44.84x10“Tg C; 耕 
地 向 林地 和 水 域 的 转换 ， 也 使 植被 WPP 有 所 增加 ; 耕地 转移 城镇 导致 植被 NPP 损失 最 大 ， 损 失 量 为 
207.76x10* Tg C， 林 地 转移 城镇 、 草 地 转 为 城镇 也 带 来 了 NPP 的 损失 。 各 类 型 土地 类 型 植被 NPP 向 城镇 
转移 ， 损 失 的 NPP 量 最 多 ， 共 损失 了 216.83x10”T g C， 因 为 城市 面积 的 不 断 增加 ， 植 被 面积 不 断 减 少 ， 
城市 的 植被 NPP 有 较 大 的 损失 。 这 表明 引起 植被 NPP 损失 的 主要 驱动 力 是 城市 化 。 

2000 一 2005 年 草地 向 耕地 转移 土地 面积 481 km*， 使 得 NPP 增加 了 272.28x10“Tg C， 而 之 后 的 五 年 
间 仅 有 3 km2 的 草地 转移 为 耕地 ， 远 远 小 于 2000 一 2005 年 转移 的 面积 , 可 能 由 于 转移 的 土地 并 未 得 到 充分 
的 耕作 ， 未 能 很 好 的 利用 ， 原 本 的 草地 荒废 ， 使 得 转移 为 耕地 的 土地 NPP 不 升 反 降 ， 导 致 了 两 个 研究 时 间 
段 相同 转移 类 型 损益 不 同 。 山 东 省 拥有 众多 的 湖泊 湿地 资源 ， 随 着 退耕 还 湿 还 湖 工程 的 实施 ， 众 多 人 工 和 
自然 湖泊 湿地 的 情况 得 到 进一步 恢复 ， 水 域 中 的 芦苇 、 薄 草 等 繁殖 能 力 强 的 植物 ， 形 成 大 范围 的 植物 群落 
( 刘 颖 ，2017) ， 良 好 的 植物 多 样 性 使 得 水 域 植被 的 净 初 级 生产 力 较 之 周边 耕地 、 草 地 的 高 ， 湖 泊 湿 地 等 
水 域 植 被 生长 状况 得 到 进一步 的 改善 。 此 外 ，2005 年 山东 省 政府 下 发 了 基本 农田 保护 工作 的 意见 ， 即 保护 
GD ”了 作物 的 生长 力 , 同时 有 效 防止 了 耕地 城镇 化 。2008 年 政府 下 发 胶东 半岛 城市 群 和 省 会 城市 群 一 体 发 展 规 
吕 划 的 通知 ， 在 经 济 发 展 的 情况 下 ， 优 先 发 展 主 城区 ， 设 立 农田 保护 区 、 水 源 保护 区 和 生态 自然 保护 区 等 ， 
一 因地制宜 ， 为 生态 系统 的 良性 发 展 提供 保障 。 


SI 4 讨论 与 结论 


一 研究 表明 ， 和 气候 因素 和 人 为 因素 等 共同 影响 着 植被 WPP 的 空间 变化 ， 面 积分 别 占 研究 区 的 9.21% 和 
一 90.79%， 说 明 人 为 因素 更 大 程度 的 影响 着 山东 省 的 植被 VPP。 对 人 类 活动 的 分 析 利 用 了 土地 覆盖 类 型 转移 

的 方法 , 其 中 城市 面积 的 增加 为 引起 植被 NPP 损失 的 主要 因素 。 文 中 非 气候 因子 驱动 型 面积 达到 90.79%， 
r= 人 类 活动 导致 的 土地 利用 类 型 转移 只 是 其 中 一 部 分 的 原因 ，2000 一 2005 年 和 2005 一 2010 年 的 土地 转移 面 
人 职 分 别 只 有 2% 和 0.92%。 人 类 活动 所 诱发 的 自然 灾害 、 大 型 工程 建设 等 因素 也 会 影响 植被 VPP 的 变化 。 
研究 区 海拔 、 坡度 的 不 同 , 使 得 气候 条 件 不 尽 相 同 , 不 同 海拔 和 坡度 所 适应 的 植被 类 型 、 土 壤 肥力 的 差别 ， 
也 会 影响 到 植被 NPP 的 大 小 ， 这 些 因素 在 今后 的 研究 中 会 进一步 讨论 分 析 。 

人 类 活动 对 植被 NPP 有 负面 减损 和 正面 增多 的 双重 作用 。 如 城市 化 进程 的 加 快 导致 对 植被 的 破坏 等 因 
素 为 负面 作用 ，2000 一 2010 年 ， 山 东 省 城市 面积 增加 了 2 787 km”， 城 市 面积 增加 导致 损失 的 NPP 共有 
323.21x10“Tg C; 而 加 大 对 农田 保护 的 投入 力度 ,加 强 农 田 水 利 基本 建设 等 措施 为 正面 影响 ， 这 十 年 间 耕 
地 面积 增加 了 788 km2， 增 加 的 耕地 NPP 共有 177.27x10“T g C， 其 中 未 利用 地 和 草地 向 耕地 的 转移 ， 以 
及 耕地 自身 生产 力 的 提高 使 得 耕地 NPP 总 量 增 加 。 

本 研究 中 ， 山 东 省 植被 NPP 整体 为 增加 趋势 ， 植 被 状态 好 转 ， 与 国内 学 者 的 研究 结果 相同 ， 苏 本 营 等 
(2010) 对 山东 省 农田 生态 系统 的 研究 发 现 ，2000 一 2006 年 农田 生态 系统 在 呈现 好 转 的 趋势 。 然 而 崔 林 丽 
等 (2016) 在 中 国 东南 部 植被 VPP 时 空格 局 变化 的 研究 中 发 现 2001 一 2010 年 东南 部 植被 NPP 整体 上 是 微 
弱 的 减少 趋势 。 由 于 植被 类 型 地形 地 貌 和 人 为 干预 等 因素 ,使 得 不 同 区 域 尺度 的 植被 NPP 存在 明显 差异 。 
研究 中 使 用 的 WPP 数据 分 辩 率 较 低 ， 无 法 精确 的 区 分 不 同 植被 类 型 , 分析 中 可 能 会 带 来 结果 的 误差 在 今 
后 的 研究 中 应 运用 更 高 分 辩 率 的 数据 ; 牛 忠 恩 等 (2016) 提 出 MOD17 产品 中 的 光 能 利用 率 导 致 其 植被 NPP 
被 低估 ， 在 今后 的 研究 中 应 结合 研究 区 实际 情况 ， 具 体 参 数 做 适当 的 调整 ， 并 结合 多 种 NPP 计算 方法 ， 提 


高 更 适用 于 区 域 研究 的 模型 精度 。 

综 上 所 述 ， 结 论 如 下 : 

(1) 2000 一 2014 年 ， 山 东 省 NPP 年 均值 为 442 gC.m”*.a1!，15 年 间 NPP 总 量 平均 为 63.16TgC， 呈 
现 先 增 加 ， 后 相对 稳定 略 有 升 高 的 趋势 。 
(2) 山东 省 植被 NPP 年 均值 呈现 明显 的 空间 分 异 特征 。 低 值 区 集中 于 济南 市 、 淄 博 市 、 莱 芜 市 、 泰 
安 市 和 济宁 市 ， 德 州 市 西南 部 、 滨 州 市 和 东营 市 北部 地 区 也 有 分 布 ， 高 值 区 位 于 东北 沿海 区 域 ， 最 大 值 在 
威海 市 。 

(3) 山东 省 NPP 增加 的 区 域 占 总 面积 的 92.36%,， 减少 区 域 占 比 为 7.64%。NPP 极 显 著 减 少 和 显著 减 
少 的 面积 合计 仅 有 1.33%, 变化 不 明显 区 域 为 62.95%, 显著 增加 区 域 19.16%，, 极 显 著 增 加 区 域 占 16.56%。 

(4) 植被 NPP 变化 的 主要 驱动 因素 是 气候 变化 和 人 类 活动 。 人 类 活动 的 影响 最 大 ， 占 研究 区 总 面积 
的 90.79%， 城 市 面积 的 增多 是 主要 影响 NPP 变化 的 因素 ， 主 要 体现 在 土地 利用 面积 中 耕地 向 城市 的 转移 
损失 了 NPP。 山 东 省 植被 NPP 受气 候 因子 的 影响 面积 占 9.21%， 其 中 降水 是 主要 的 影响 因素 。 


~ 
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